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New astrophysical instruments such as skA (square kilometer Array) 
and IXO (formerly Constellation X) promise the discovery of tens of thou- 
sands of new isolated rotating neutron stars (pulsars), neutron stars in 
low- mass X-ray binaries (LMXBs), anomalous X-ray pulsars (AXPs), and 
soft gamma repeaters (SGRs). Many of these neutron stars will experi- 
ence dramatic density changes over their active lifetimes, driven by either 
stellar spin-up or spin-down, which may trigger phase transitions in their 
dense baryonic cores. More than that, accretion of matter onto neutron 
stars in LMXBs is believed to cause pycno-nuclear fusion reactions in the 
inner crusts of neutron stars. The associated reaction rates may be dras- 
tically altered if strange quark matter would be absolutely stable. This 
paper outlines the investigative steps that need to be performed in order to 
explore the thermal response of neutron stars to rotationally-driven phase 
transitions in their cores as well as to nuclear burning scenarios in their 
crusts. Such research complements the exploration of the phase diagram of 
dense baryonic matter through particle collider experiments, as performed 
at RHIC in the USA and as planned at the future Facility for Antiproton 
and Ion Research (FAIR) in Darmstadt, Germany. 

PACS numbers: 26.60+c, 97.60.Gb, 97.60.Jd 

1. Introduction 

On the Earth, particle collider experiments enable physicists to cast a 
brief glance at the properties of ultra-dense and hot baryonic matter (Fig. 
[1]). On the other hand, it is estimated that galaxies like our Milky Way 
contain between 10 8 and 10 10 collapsed stars known as neutron stars, which 
harbor ultra-dense matter permanently in their cores. This key feature to- 
gether with the unprecedented progress in observational astrophysics, which 
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Fig. 1. Schematic phase diagram of strongly interacting matter [17]. The dashed 
circle indicates the portion of the phase diagram that can be explored by studying 
either hot and newly formed proto neutron stars (arrow pointing downwards), 
or cold rotating neutron stars (vertical double-headed arrow) such as millisecond 
pulsars and neutron stars in low-mass X-ray binaries (LMXBs). 



is expected to be excelled by future observatories such as the square kilo- 
meter Array (skA) and Constellation-X, make neutron stars superb astro- 
physical laboratories for a wide range of physical studies. These studies 
concern nuclear fusion processes on the stellar surface, pycnonuclear reac- 
tions in electron degenerate matter at sub-nuclear densities, and the possible 
formation of boson condensates and other novel states of baryonic matter- 
like color superconducting quark matter-at super-nuclear densities. (For 
overviews, see, for instance [HOEIIHEIEIEIEIEI HO]-) More than that, 
there is the very intriguing theoretical suggestion that strange quark matter 
could be more stable than atomic nuclei, known as the strange quark matter 
hypothesis, in which case neutron stars should be largely composed of pure 
strange quark matter [T2j [13] . If quark matter exists in neutron stars it 
ought to be a color superconductor [5] [HI [151 [16] . Other testable implica- 
tions of the strange quark matter hypothesis concern the possible existence 
of a new class of white dwarfs, known as strange white dwarfs |18[ 119]. and 
the drastic alteration of heavy-ion reaction rates in the deep crustal layers 
of neutron stars [20], if strange quark matter nuggets should be present in 
these layers. 



2. Phase transitions in the cores of neutron stars 



Rotating neutron stars are called pulsars. Three distinct classes of pul- 
sars are currently known. These are (1) rotation powered pulsars, where the 
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Fig. 2. Composition of a rotating neutron star in equatorial direction (left panel) 
and polar direction (right panel) |21j . The star's mass at zero rotation is 1.40 Mg. 



loss of rotational energy powers the emitted electromagnetic radiation, (2) 
accretion-powered (X-ray) pulsars, where the gravitational potential energy 
of the matter accreted from a low-mass companion is the energy source, 
and (3) magnetars, where the decay of a ultra-strong magnetic field pow- 
ers the radiation. Depending on star mass and rotational frequency, the 
matter in the core regions of neutron stars may be compressed to densities 
that are up to an order of magnitude greater than the density of ordinary 
atomic nuclei. This extreme compression provides a high-pressure envi- 
ronment in which numerous subatomic particle processes are likely to take 
place [TJ [H[3l [9]. The most spectacular ones stretch from the generation of 
hyperons and baryon resonances (E, A, E, A), to quark (u,d,s) deconfine- 
ment, to the formation of boson condensates (ir~ , K~ , H) [H [2j HI [61 [8] [9] . 
Rapid rotation can change the structure and composition of neutron stars 
substantially, depending on the equation of state of ultra-dense baryonic 
matter. The most rapidly rotating, currently known neutron star is pulsar 
PSR J1748-2446ad, which rotates at a period of 1.39 ms (which corresponds 
to a rotational frequency of 719 Hz) [25]. It is followed by PSRs B1937+21 
[23j and B1957+20 [SJj whose rotational periods are 1.58 ms (633 Hz) and 
1.61 ms (621 Hz), respectively. Finally, we mention the recent discovery of 
X-ray burst oscillations from the neutron star X-ray transient XTE J 1739- 
285 [251 Eg], wnich could suggest that XTE J1739-285 contains an ultrafast 
neutron star rotating 1122 Hz. Rotating neutron stars appear as much 
better probes for the structure of dense baryonic matter than non-rotating 
neutron stars, primarily because the particle compositions in rotating neu- 
tron stars are not frozen in, as it is the case for non-rotating neutron stars, 
but are varying with time. The associated density changes can be as large 
60% [2J in neutron stars in binary stellar systems (e.g., LMXBs), which are 
being spun up to high rotational frequencies, or isolated rotating neutron 
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Fig. 3. Same as Figure [H but for a non-rotating stellar mass of 1.70 M© |27j . 



stars (e.g. isolated millisecond pulsars) which are spinning down to low fre- 
quencies because of the emission of gravitational radiation, electromagnetic 
dipole radiation, and a wind of electron-positron pairs. As an example, we 
show the rotationally-driven restructuring effects in the cores of standard 
neutron stars in Figs. [2] and [3l Qualitatively similar restructuring effects 
were obtained for neutron stars containing quark matter [2j |6] . 



3. Thermal evolution of rotating neutron stars 

Figure H] illustrates schematically the differences of the cooling behav- 
ior of rotating and non-rotating neutron stars. Computing the thermal 
evolution of rotating neutron stars, however, is considerably more compli- 
cated than computing the cooling behavior of non-rotating neutron stars, 
for several reasons [2]. First, stellar rotation requires solving Einstein's 
field equations for rotationally deformed fluid distributions, which renders 
the problem 2-dimensional. Second, the general relativistic frame dragging 
(Lense-Thirring) leads to the appearance of additional non-linear field equa- 
tion. Third, the determination of the general relativistic Kepler frequency, 
which sets an absolute limit on stable rapid rotation, leads to an additional 
self-consistency condition [2j. Fourth, the thermal transport equations need 
to be solved for general relativistic, non-spherical fluids that may experience 
anisotropic heat transport. Because of all these complications, a fully self- 
consistent general relativistic treatment of the cooling of rotating neutron 
stars has not been achieved yet. However, first steps toward this goal were 
made in [28]. The basic cooling features of a neutron star are easily grasped 
by considering the energy conservation relation of the star in the Newtonian 
limit [22]. This equation is given by 



dE th /dt = C v dT/dt = —L u -L^ + H, 



(1) 
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Fig. 4. Schematic difference of the cooling behavior of rotating and non-rotating 
neutron stars. The composition is frozen in in non-rotating neutron stars, but 
varies strongly with rotational frequency, v, in rotating neutron stars during spin- 
down/up. This restructuring effect can alter the cooling behavior of rotating neu- 
tron stars significantly. Further differences arise from the different geometry (met- 
ric) of rotating stars and from anisotropic heat transport. 



where is the thermal energy content of a neutron star, T its internal 
temperature, and Cy its total specific heat. The energy sinks are the total 
neutrino luminosity, L u , and the surface photon luminosity, L 7 . The source 
term H includes all possible heating mechanisms [29], which, for instance, 
convert magnetic or rotational energy into heat. The dominant contribu- 
tions to Cy come from the core whose constituents are leptons, baryons, 
boson condensates and possibly deconfined superconducting quarks. When 
baryons and quarks become paired, their contribution to Cy is strongly sup- 
pressed at temperatures smaller than the critical temperatures associated 
with these pairing phases. The crustal contribution is in principle domi- 
nated by the free neutrons in the inner stellar crust but, since these are 
extensively paired, practically only the nuclear lattice and electrons con- 
tribute. Extensive baryon, and quark, pairing can thus significantly reduce 
Cy. In order to derive the general relativistic version of Eq. ([1]) for rotat- 
ing stars, one needs to solve Einstein's field equations using the metric of a 
rotationally deformed fluid [2], 



where e 2 ^ = e 2(a+/3 V 2 sin 2 8 and e 2u} = e 2(a . The quantities i/, 4> and 
u denote metric functions, and N v accounts for frame dragging caused by 
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the rotating fluid. All these functions are to be computed self-consistently 
from Einstein's field equation, G a @ = 87rT a ^, where T a @ denotes the fluid's 
energy momentum tensor. In the case of neutron star cooling, the equations 
which describe the conservation of energy and momentum in a co-moving 
reference frame follow from T a P.g = q a , where 



nCt/3 



(J H r Hq 
H r J/3 
H e J/3 












J/3 / 



-H r /A 
-H e /A 
V -H v /A ) 



(3) 



The quantity q a describes the interaction of photons with matter, J = aT 4 
is the energy density of the photon field, H$ is the i th component of the heat 
flux, e the emitted energy, and A the mean free path of photons. From Eq. 
([3]) one obtains 



q a = e P_ T& p^ +l a_ r fP + 1 P_ T aj i 
where the 7's denote the Ricci rotation coefficients given by 7?. 



(4) 



Combining Eqs. ([3]) and ^ leads to the general relativistic transport equa- 
tions 
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(8) 



where H { = re 2u+<t>+u) Hi/Y , f = e u T/Y, and the Lorentz factor T = (1 - 
U 2 )^ 1 / 2 . In the case of uniform rotation, £1 =const, Eqs. © to ([8]) reduce 
to 
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The standard cooling equations of spherically symmetric, non-rotating neu- 
tron stars are obtained from Eqs. ([9]) to (fTTj) for Q = and dgT = [2]. 
This project aims at solving Eqs. to dTJ) for the temperature distri- 
bution T(r, 9; t) of non-spherical, (possibly differentially) rotating neutron 
stars. The boundary condition are given by defining Hf at r = 0, R, and 
Hg at 9 = 0, 7r/2 and at r = R, with R denoting the stellar radius. The 
choice for the star's initial temperature, T(r, 9; t = 0), is typically chosen as 
f = 10 11 K. Equations © and © can be solved for Hf and Hg, differen- 
tiated with respect to d r and dg, respectively, and substituted into Eq. ([5]), 
which leads to the following parabolic differential equation, 

8 t f = -± e 2,/ -^-r S m9Ue u+ ~<-Z-^(d r ri + -d e n) 
r Cy Cy V r J 

+ 1 n -e 3u ~^ -j-fd r f r 2 Ksm9e" / ( d r f + T 2 Ue- 2u+ ^fd r Q)) 
r z sin 9 T Cy V V v 

+ (r 2 Ksin0e 7 (a e r + r 2 C/ e -^ +7 Ta^))) , (12) 

with the definitions rsin^e^^^ 7 = and e _l/+ ^ = e a ~@. It is this dif- 
ferential equation that needs be solved numerically in combination with a 
general relativistic stellar rotation code, as illustrated schematically in Fig. 
[5j A 2-dimensional, general relativistic stellar rotation code is to be used 
to determine the metric functions, frame dragging frequency, pressure and 
density gradients, and particle compositions of a deformed neutron star as 
a function of rotational frequency. Depending on the cooling channels that 
are active at a given rotational frequency, the numerical outcome obtained 
from the rotation code then serves as an input for the general relativistic 
2-dimensional thermal evolution code, which is used to determine the lumi- 
nosity, and thus the surface temperature, of a deformed rotating neutron 
star. As outlined in Sect. [21 because of stellar spin-down or spin-up, the 
density in such stars may change dramatically so that new cooling chan- 
nels open up (or close) with time (Fig. SJ). The new cooling channels imply 
neutrino emissivities, heat capacities, and thermal conductivities that are 
different from the original ones, thus altering the stars thermal response. It 
is this response in the thermal behavior of rotating neutron stars that car- 
ries information about the properties of the matter in the dense baryonic 
stellar cores and the deep crustal layers. Several key research activities that 
can be addressed with these codes are briefly summarized next: 



Study the impact of latent heat, released by rotation-driven quark- 
hadron phase transitions in the cores of isolated rotating neutron stars, 
on the thermal evolution of such objects |28j. 
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, 5. Combination of numerical codes required to perform thermal evolution 
ies of rotating neutron stars. 



• Simulate the impact of latent heat release on the thermal evolution 
of accreting neutron stars in LMXBs. Because of the spin-up of such 
neutron stars, these objects becomes progressively more decompressed 
with time, which may lead to drastic changes in their core composi- 
tions, causing the destruction of certain states (K— condensate, quark 
matter) of dense baryonic matter that were initially present at slow 
stellar rotation rates. This modifies the star's thermal evolution and 
may lead to signals indicative of the existence of phase transitions in 
the star. 

• Study the impact of energy released by pycnonuclear burning of mat- 
ter in the inner crusts of neutron stars. In particular, the dependence 
of the stellar surface temperature on the possible presence of strange 
quark matter nuggets in the inner crusts of neutron stars can be ex- 
plored. As shown in Fig. [61 the nuclear reactions rates are strongly 
modified if strange quark matter nuggets would be present in the 
crusts of neutron stars [2D] . 

• The thermal evolution of neutron stars undergoing episodes of intense 
accretion, alternated by long periods of quiescence, i.e. soft X-ray 
transients (SXRTs) such as SAX J1808. 4-3658, can be explored. This 
neutron star is the first known accretion powered millisecond pulsars 
rotating rapidly (at 2.5 ms, 400 Hz). At such small spin periods 
certain phases of dense matter may have been spun out of the core, 



paper'FWeber printed on November 26, 2009 



9 




altering the neutrino loss. The quiescent emission of SAX J1808.4- 
3658, therefore, opens the possibility of using this object as a tool 
for probing the core compositions of neutron stars |29(, 130] . Another, 
and still poorly explored, aspect of SXRTs is their short time-scale 
thermal response to the accretion phases. This is an aspect of the 
problem from which extremely important information about both the 
structure of the neutron star crust and the thermal state of its core 
can be obtained |314 I32j . and about which intriguing observational 
results were found [33]. 

• Another interesting problem concerns the thermal evolution of neu- 
tron stars immediately after they went through the proto-neutron star 
phase. Of particular interest is the crystallization behavior of such 
matter and the occurrence of thermoelectric instabilities at the core- 
crust boundary. 
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